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Wc present the first statistically meaningful results from two-K^ interfcromctry in heavy-ion 
collisions. A model that takes the effect of the strong interaction into account has been used to fit 
the measured correlation function. The effects of single and coupled channel were explored. At the 
mean transverse mass (mr) = 1.07 GeV, we obtain the values R = 4.09 ± 0A6{stat.) ± 0.31{sys) 
fm and A = 0.92 ± 0.23{stat) ± 0.13{sys), where R and A are the invariant radius and chaoticity 
parameters respectively. The results are qualitatively consistent with mr systematics established 
with pions in a scenario characterized by a strong collective flow. 
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I. INTRODUCTION 

Lattice QCD calculations predict that a phase tran- 
sition from hadronic matter to a new state of matter 
called a Quark Gluon Plasma (QGP) occurs at suffi- 
ciently large energy densities [ij. Creation and study 
of such a de-confined state of matter is the primary goal 
of the heavy-ion collisions program at the Relativistic 
Heavy-Ion Collider (RHIC). A first order phase transi- 
tion from the QGP back to normal hadronic matter is 
believed to delay the expansion of the hot reaction zone 
created in the collision [J. A delayed expansion means 
a long duration of particle emission, leading to a large 
source size. 

The measurement of the space-time extent of the par- 
ticle emitting region has been one of the important goals 
in high energy experiments for several decades |E |j, |^ . 
These measurements are based on the sensitivity of parti- 
cle momentum correlations to the space-time separation 
of the particle emitters due to the effects of quantum 
statistics (QS) and final state interaction(FSI). For iden- 
tical particles, the QS symmetrization (antisymmetriza- 
tion) is usually the dominant source of the correlation 
and, due to the interference of the amplitudes corre- 
sponding to various permutations of identical particles, 
this measurement is often called particle interferometry 
(See reference for a review). 

Most of the particles produced in relativistic heavy- ion 
collisions are pions and, as a result, pion interferometry 
has been a particularly useful tool in correlation studies. 
High statistics data from colliders like RHIC have also 
made it possible to study kaon correlations. In this Let- 
ter, we present the first results on two-K^ correlations in 
central Au-Au collisions at y^sWW = 200 GeV measured 
by the STAR(Solenoidal Tracker at RHIC) experiment 
at RHIC . 

It is known that a significant fraction of pions come 
from resonance decays after freeze-out thus complicating 
the pion interferometry measurements. While the direct 
pion source may be inherently non-Gaussian, the reso- 
nances extend the source size due to their finite lifetime, 
introduce an additional essentially non-Gaussian distor- 
tion in the two-pion correlator and reduce the fitted cor- 
relation strength. Due to the limited decay momenta, the 
decay pions populate mainly the low momentum region, 
thus introducing an additional pair momentum depen- 
dence in the correlator. 

Kaon interferometry, on the other hand, suffers less 
from resonance contributions and could provide a cleaner 
signal for correlation studies than pions 0, II]- Higher 
multi-particle correlation effects, that might play a role 
for pions, should be of minor importance for kaons since 
the kaon density is considerably smaller than the pion 
density at RHIC (^/i^r^r = 200 GeV). The pion mul- 
tiplicity has increased by approximately 70% from the 
SPS(Vs^^ = 17.3 GeV) to RHIC 9]. The interferome- 
try radii however remain almost the same [T(l ITll | . The 
strangeness distillation mechanism might further in- 



crease any time delay QGP signature. This mechanism 
could lead to strong temporal emission asymmetries be- 
tween kaons and anti-kaons |l3j| . thus probing the latent 
heat of the phase transition. 

Particle identification for pions, via their specific ion- 
ization (energy loss per unit length or dE/dx), works only 
up to about 700 MeV/c. Neutral kaons, on the other 
hand, can be identified up to much higher momentum 
using their decay topology. This allows for the extension 
of the interferometry systematics to a higher momentum 
than is presently achievable with pions, and thus provide 
a means to probe the earlier times of the collision. The 
effect of two-track resolution, which is a limiting factor 
in charged particle correlations, is also small. The ab- 
sence of Coulomb FSI suppression together with small 
contributions from resonance decays make neutral kaon 
correlations a powerful tool to investigate the space time 
structure of the particle emitting source. 

The OPAL J^J and ALEPH collaborations have 
measured correlations of neutral kaons from hadronic de- 
cays of Z° in e+e" collisions at LEP. The WA97 experi- 
ment at CERN (1^] attempted to measure K^K^ correla- 
tions but did not see a significant enhancement of neutral 
kaon pairs in the region of small momentum difference 
due to a lack of sufficient statistics. 



II. THE STAR EXPERIMENT 

The STAR detector consists of several detector 
subsystems in a large solenoidal magnet that provides 
a uniform 0.5 Tesla field. For the data used in this 
analysis, the main s etup consisted of the time projec- 
tion chamber (TPC) [ig for charged particle tracking, a 
scintillator trigger barrel (CTB) surrounding the TPC for 
measuring charged particle multiplicity, and two zero de- 
gree calorimeters (ZDC) 19] located upstream and down- 
stream along the axis of the TPC and beams to detect 
spectator neutrons. With full azimuthal coverage over 
]?7] < 1 and an almost 100% efficiency for minimum ion- 
izing particles, the CTB provides a good estimate of the 
number of charged particles produced in the mid-rapidity 
region. The number of neutrons detected in the ZDC's is 
identified with the amount of energy deposited in them. 
The collision centrality is determined by correlating the 
energy deposition in the ZDC with the number of mini- 
mum ionizing particles detected by the CTB. 



A. Data Selection 

For this analysis, events from the ZDC central trig- 
ger (0 — 10% of the total hadronic cross section) were 
used with an event vertex within ±25 cm of the center of 
the TPC, along the beam axis. Approximately 2.5 x 10^ 
events with about 3 K'^ per event on the average were 
analyzed. Here we discuss /Cg-specific issues only, as de- 
tails of pion interferometry at the STAR experiment have 



4 



o 

1 140 



= 100 
III 



xlO 



80 



60 



40 



20 



: t 


\ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


\ 



0.4 0.42 0.44 0.46 



0.5 0.52 0.54 0.56 0.58 
Mass(GeV/c') 



FIG. 1: The K° invariant mass distribution. The range in 
transverse momentum is from 0.5 GeV/c to 3.5 GeV/c and 
rapidity is between -1.5 and 1.5. Kaon candidates falhng in 
the mass ranee from 0.48 GeV/c^ to 0.51 GeV/c^ indicated 
by the shaded region, were selected for this correlation study. 
The corresponding mass is 495.6 ±6.8 MeV/c^. 



been discussed in |23]. The has a mean decay length 
(ct) of 2.7 cm and decays via the weak interaction into 
TT"*" and TT^ with a branching ratio of about 68%. The 
mass and kinematic properties of the are determined 
from the decay vertex geometry and daughter particle 
kinematics |21|. Neutral kaon candidates are formed out 
of a pair of positive and negative tracks whose trajecto- 
ries point to a common secondary decay vertex which is 
well separated from the primary event vertex. All neutral 
kaon candidates, with invariant masses from 0.48 GeV/c^ 
to 0.51 GeV/c^, transverse momentum from 0.5 GeV/c 
to 3.5 GeV/c and rapidity between -1.5 and 1.5 have been 
considered. The daughter particle tracks are required to 
have a minimum of 15 TPC hits and a distance of closest 
approach to the primary vertex greater than 1.3 cm. 



B. The Correlation function 

Experimentally, the two-particle correlation function is 
defined as 



C2(Q) 



AjQ) 
BiQY 



(1) 



where A{Q) represents the distribution of the invariant 
relative momentum Q — \/~q^q^n = Pi — P21 for a 
pair of particles from the same event. The possibility of 
a single neutral kaon being correlated with itself, i.e., cor- 
relation between a real if ^ and a fake reconstructed 
from a pair which shares a daughter of the real K^, was 
eliminated by requiring that kaons in a pair have unique 



daughters. We have also explored effects from splitting 
of daughter tracks by looking at the angular correlation 
between the normal vectors to the decay planes of the 
in a given pair. No enhancement at very small angles 
was observed indicating no significant problem from track 
splitting. B{Q) is the reference distribution constructed 
by mixing particles from different events with similar Z- 
vertex positions (relative z position within 5 cm). The 
individual for a given mixed pair are required to pass 
the same single particle cuts applied to those that go into 
the real pairs. The mixed pairs are also required to sat- 
isfy the same pairwise cuts applied to the real pairs from 
one event. The efficiency and acceptance effects cancel 
out in the ratio ^^gj . 



C. Data Analysis 

Figure ^ shows the invariant mass distribution of the 
neutral kaons based on the set of cuts described above. 
The background is characterized by a polynomial fit to 
the distribution outside the mass peak. The observed 
mass 495.6 ± 6.8 MeV/c^ is consistent with the accepted 
value ■ The signal and background for the mass range 
from 0.48 GeV/c^ to 0.51 GeV/c^ considered in this 
analysis are shown by the shaded regions. 
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FIG. 2: The signal to (signal-|-background) ratio as a 
function of the transverse momentum pt- The data points 
correspond to a decay length (DL) greater than 6 cm. The 
kaons selected fall in the mass range from 0.48 GeV/c^ to 
0.51 GeV/c^ which is also the mass range for the correlation 
analysis. The errors are only statistcal. 

After tuning several kinematical and detector related 
cuts to remove most of the background, some residual 
noise still remains. This calls for a knowledge of the 
signal to background ratio within the selected invariant 
mass range to make corrections to the measured correla- 
tion function. For neutral kaons, the decay length (DL) 
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FIG. 3: The Kt distribution of the K° pairs. The range 
in transverse momentum of the single particles is from 0.5 
GeV/c to 3.5 GeV/c. The distribution in (a) corresponds to 
Q < 0.2 GeV/c and that in (b) is for Q < 0.1 GeV/c, i.e., (b) 
is a subset of (a). The two histograms in each panel are for 
low (dashed) and high (full) pair purity.' 



and distance of closest approach (DCA) to the interaction 
vertex were two of the parameters for which it was diffi- 
cult to determine where to apply the cuts. Various DCA 
and DL cut combinations were investigated by varying 
the DCA from 0.3 cm to 0.8 cm in steps of 0.1 cm and 
the DL from 2.0 cm to 6.0 cm in steps of 1.0 cm. Figure 
12 displays an example of the signal to background ratio 
as a function of pr for DL > 6 cm and various DCA val- 
ues. The single particle purity gets worse as the DCA 
gets larger for the given DL cut. If one instead looks at 
a fixed DCA and varies the DL cut instead, the purity 
gets better with increasing decay length. 

The effect of momentum resolution on the correla- 
tion functions has also been investigated using simulated 
tracks from K'^ decays with known momenta, pi„, embed- 
ded into real events. The reconstructed momenta of the 
embedded tracks, pVec, are then compared with pin- The 
distributions of with respect to Pm are then fit 

to Gaussians to obtain the RMS widths. These are used 
to characterize the momentum resolution of the detector. 
The resolution in p lies between 1% and 2% for the pt 
range used in this analysis. 



The top panel in Figure |3| shows the Kt distribution 
for Q < 0.2 GeV/c where Kt = {\pit + P2t|)/2. The 
correponding number of pairs for the distribution with 
low pair purity is approximately 1.92 x lO'* and that for 
the one with the high pair purity is about 5.5 x 10"^. 
The distribution in the bottom panel corresponds to pairs 
with Q < 0.1 GeV/c, with 2.7 x 10^ pairs for the low 
pair purity distribution and 7.8 x 10^ for the high pair 
purity distribution. It is clear that the shape of the Kt 
distribution changes with the pair purity and, as a result, 
so does (Kt), the mean of the distribution. The mean 
Kt varies almost linearly with pair purity. For the lowest 
pair purity value of « 52%, (Kt) « 0.805 GeV/c. At the 
highest pair purity value of « 89%, (Kt) « 1.07 GeV/c. 
The dependence of (Kt) on the pair purity together with 
the fact that the radii may vary with Kt implies that 
varying the pair purity may change the measured radii. 
In this analysis, the correlation function is integrated over 
all Kt since the statistics are not sufficient to make a Kt 
dependent study. 

Corrections to the raw correlation functions were ap- 
plied according to the expression 

„ / „N _ CmeasurediQ) ~ 1 /„n 

'-^correctedK'^ j ~ r, ■ D -j- 1 ' ^ ' 

FairFurity[Li) 

where the pair purity was calculated as the product of 
the signal(S') to signal plus background [S + B) ratios of 
the two K^ of the pair (i,j) 

^ -Xpu) ^^{Pt,) (3) 



PairPurity{Q) 



S + B"^"' S + B' 

The pair purity, P air Purity [Q), has been found to 
be independent of Q over the range of invariant four- 
momentum difference considered. As a result, an average 
value over Q of the pair purity has been used to correct 
the correlation function for each set of cuts considered. 

Figure 0] shows the experimental K^K^ correlation 
function before and after corrections for purity and mo- 
mentum resolution are applied. It can be seen that the 
effect of momentum resolution is comparable to that of 
purity correction. The one-dimensional correlation func- 
tion is usually fitted to a Gaussian 



C(Q) = • (1 + A • e 



) 



(4) 



where N and R are respectively the normalization and 
size parameter, the latter characterizing the width of the 
Gaussian distribution of the vector r* of the relative dis- 
tance between particle emission points in the pair c.m.s.: 

(5) 



(Pi 



oc e 



The parameter A measures the correlation strength. In 
the absence of FSI, A equals unity for a fully chaotic 
Gaussian source, up to a suppression due to the kaon 
impurity and finite momentum resolution. Theoretically, 
it can be less than unity due to partial coherence of the 
kaon field, resonance decays and the non-Gaussian form 
of the correlation function. Also neglecting FSI can affect 
(suppress or enhance) the value of this parameter. 
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FIG. 4: The K^K° correlation function. The solid circles 
are for uncorrected data. The squares correspond to the case 
where the data have been corrected for pair purity. The tri- 
angles represent the data after correcting for pair purity and 
momentum resolution. The errors are ststistical only. 



III. FINAL STATE INTERACTION IN THE 
NEUTRAL KAON SYSTEM 

The production of the neutral kaon system, and 
, is attributed to the strong interaction which con- 
serves the strangeness quantum number. An interesting 
property of neutral kaons is that the can change into 
a if through a second order weak interaction. However, 
the particles that we normally observe through weak de- 
cay channels in the laboratory are not K'^ and K'^ |23l |. 
Neglecting the effects of CP violation, the observed weak 
interaction eigenstates are given by 



1 

1 

72 



(6) 



where \K^) and l-K^f) are the state vectors of the short 
and long lived neutral kaons, to which experiments have 
access via measurements of their decay products, which 
are mainly pions. The state vector of the K^K^ system 
is then given by the expression 



(7) 



Now, if a K'^K° pair comes from K°K° K°), it is 
subject to Bose-Einstein (BE) enhancement as it origi- 
nates from an identical boson pair. On the other hand, 
the K'^ and K'^ are two different particles and one may 



not expect correlations if one comes from and 
the other one from K^. Nevertheless, it can be shown 
|23 | (see also jHllEllS) that only the symmetric part 
of the K^K^ amplitude contributes to the K^K'^ sys- 
tem and thus also leads to a Bose-Einstein enhancement 
at small relative momentum (on the contrary, only the 
anti-symmetric part of the K^K^ amplitude contributes 
to the K^Ki system and leads to the "Fermi-Dirac like" 
suppression). The K^K^ correlation thus includes a 
unique interference term that may provide additional 
space-time information. Here only the K^K^ correla- 
tion is considered since most of the if ^ decay outside the 
STAR TPC and are not accessible. 

The strong FSI has an important effect on neutral kaon 
correlations due to the near threshold resonances, /o(980) 
and_ao(980) These resonances contribute to the 

K^K'^ channel and lead to the s-wave scattering length 
dominated by the imaginary part of ^1 fm. Based on the 
predictions of chiral perturbation theory for pions p9| 
the non-resonant s-wave scattering lengths are expected 
to be ~0.1 fm for both K^K^ and K^K^ channels and 
can be neglected to a first approximation. 

To calculate the K^K^ correlation function, we as- 
sume isT^'s and K^^s emitted by independent single-kaon 
sources so that the fraction of K^K^ pairs originating 
from K^K^ system is a = (1 — e^)/2, where e is the 
K^-K'^ abundance asymmetry. We have put a = 1/2 
based on the negligible K^-K^ abundance asymmetry 
of 0.018 ± 0.106 as measured under the same conditions 
by the STAR experiment (s^. The correlation func- 
tion is calculated as a mixture of the average squares 
of the properly symmetrized K^K^, K^K^ and non- 
symmetrized K^K^ wave functions, weighted by the re- 
spective K^K^ fractions. To average over the relative 
distance vector^, we use the Lednicky & Lyuboshitz an- 
alytical model [23, assuming r* is distributed according 
to Eq. jSJ with a Gaussian radius R. The model assumes 
that the non-symmetrized wave functions '^_j:,ir*) de- 
scribing the elastic transitions can be written as a su- 
perposition of the plane and spherical waves, the latter 
being dominated by the s-wave, 



fik* 



Ak* r* 



(8) 



where k* = Q/2 is the three-momentum of one of the 
kaons in the pair rest frame and f{k*) is the s-wave scat- 
tering amplitude for a given syst_em^ Neglecting the scat- 
tered waves for the K^K'^ and K'^K^ systems (the corre- 
sponding f{k*) = 0) one obtains the following expression 
for the K2k2 correlation function 



C{Q) = l + e-'^"^ 
mj{k*) 



fik*) 



R 



FiiQR)~^^^F,{QR) 



, (9) 



0Fi? R 
where Fi{z) = dxe^ /z and F2{z) = (1 — e^^ )/z 
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The s-wave K^K'^ scattering amplitude /(fc*) is domi- 
nated by the near threshold s-wave isoscalar and isovec- 
tor resonances /o(980) and ao(980) characterized by their 
masses rrir and respective couplings 7^ and 7'. to the 
KKjTTTT and KK,ttt] channels. Associating the ampli- 
tudes // at isospin / = and I — I with the resonances 
f = fo a-nd ao respectively, one can write |28l l32| 



/(fc*)-[/o(fc*) + /i(fc*)]/2, 



(10) 



(11) 



Here s = 4(m|- -I- fc*^) and k'^ denotes the momentum 
in the second (tttt or ttt]) channel with the corresponding 
partial width FJ, — j'^k'^/rrir. 

There is a great deal of uncertainty in the properties 
of these resonances due to insufficiently accurate exper- 
imental data and the different approaches used in their 
analysis. Fortunately, the dominant imaginary part of 
the scattering amplitude is basically determined by the 
ratios of the JqKK to /otttt &u.d anKK to ao-Krj couplings 
whose variation is rather small |33l | . In this paper we use 
the resonance masses and couplings from (a) Martin et 
al. m, (b) Antonelh [2I, (c) Achasov et al. (d) 
Achasov et al. [1^ (see Table P) to demonstrate the im- 
pact of their characteristic uncertainties on the calculated 
correlation function. 



Ref. 








rUao 


lanKK 


laoTTV 


a 


0.978 


0.792 


0.199 


0.974 


0.333 


0.222 


b 


0.973 


2.763 


0.5283 


0.985 


0.4038 


0.3711 


c 


0.996 


1.305 


0.2684 


0.992 


0.5555 


0.4401 


d 


0.996 


1.305 


0.2684 


1.003 


0.8365 


0.4580 



and i? = 3 fm as input radii with the normalization factor 
N and A both set to unity. 



Mass and coupling parameters 
from row (a) of Table I 

R = 6 fm 

R = 3fm 

Mass and coupling parameters 
from row (d) of Table I 

R = 6 fm 

R = 3fm 




0.02 0.04 0.06 0.08 



0.12 0.14 0.16 0.18 0.2 
Q (GeV/c) 



FIG. 5: Theoretical correlation functions for input Gaussian 
sources of i? = 6 fm and 7? = 3 fm with A = 1, A'' = 1 The 
resonance masses and coupling constants are from Table Q 

The results indicate that the effect of the strong FSI 
in the K^K'^ system is to give rise to a repulsive-like 
component causing the correlation function to go below 
unity. 



IV. EXPERIMENTAL RESULTS 



TABLE I: The /o and ao masses and coupling parameters, all 
in GeV, from (a) Martin et al. (b) Antonelh et a/.Q, 

(c) Achasov et aZ.js^ and (d) Achasov et al. j3^. 

We have taken into account the normalization and cor- 
relation strength parameters N and A by the substitution 
C{Q) ^ N -[X- C{Q) + (1 - A)] . Following Ref. Hi, 
we have also included a small contribution of the inelas- 
tic transition between the coupled channels K^K~{= 2) 
and K^K^{= 1) (see Appendix for more details). Be- 
sides a direct contribution of the average square of the 
corresponding wave function ^'^^^j. {?*) given in Eq. (|12|l . 
this transition also leads to a modification of the ampli- 
tude f{k*) in the wave function of the elastic transition 
in Eq. (|SJ). Instead of Eq. (|1U|I . this amplitude is now rep- 
resented by the element /^^ of a 2 x 2 matrix /c defined 
in Eq. H13I) . We have further considered the correction 
ACj^jf in Eq. (|16|1 due to the deviation of the spherical 
waves from the true scattered waves in the inner region 
of the short-range potential, which is of comparable size 
to the effect of the second channel. 

Figure |31 shows the theoretical correlation functions for 
two sets of resonance parameters from TableQ]with R = 6 



The experimental correlation functions are fit using the 
Lednicky & Lyuboshitz |2^] model to take into account 
the effect of the strong FSI. The free parameters are the 
radius R characterizing the separation r* of the particle 
emission points in the pair rest frame, the normalization 
N, and A. This fitting was done assuming the Gaussian 
r* -distribution of Eq. (|SJ). 

The fit results are summarized in Table ^ for var- 
ious sets of resonance parameters. The normalization 
N — 1.03 in all cases. The difference between the sin- 
gle channel and coupled channel fits is very small, but 
it is the coupled channel fit results which are more ac- 
curate. Figure shows an example of the model fits to 
the experimental correlation function. A Gaussian fit to 
the correlation function gives R = 5.02 ± 0.61 fm and 
A = 1.08 ± 0.29. One can see that a Gaussian fit can- 
not account for the C{Q) < 1 part of the data which 
are fit better if the strong FSI is included. Figures [7| 
and IHl show the dependence of the extracted Rinv and A 
parameters as a function of the Pair Purity before, (a), 
and after, (b), correcting for this impurity. The data 
points are not independent of each other as a low purity 
data may contain some or all of the high purity data. 
The fit results are not sensitive to the resonance param- 
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3 
U 

2.5 



Gaussian fit to data 

Fit to theory with mass and 
coupling constants from Table I 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

Q (GeV/c) 

FIG. 6: Fits to experimental correlation function including 
the strong interaction with resonance masses and coupling 
constants from Table |I] The corresponding /DOF are (a) 
1.053, (b) 1.048, (c) 1.045 and (d) 1.046. A simple Gaussian 
fit, with / DOF = 0.816, is also shown for comparison. The 
errors are only statistical. 



Rinv (fm) 


1-ch. fit 


2-ch. fit 


a 


3.90 ± 0.45 ± 0.37 


4.07 ±0.46 ±0.31 


b 


3.89 ± 0.44 ± 0.35 


4.09 ±0.46 ±0.31 


c 


3.96 ± 0.45 ± 0.34 


4.14 ±0.47 ±0.31 


d 


3.91 ±0.44 ±0.34 


4.07 ±0.45 ±0.29 


X 


1-ch. fit 


2-ch. fit 


a 


0.89 ±0.21 ±0.10 


0.98 ±0.24 ±0.14 


b 


0.83 ±0.20 ±0.10 


0.93 ±0.23 ±0.13 


c 


0.81 ±0.20 ±0.09 


0.90 ±0.23 ±0.12 


d 


0.78 ±0.19 ±0.09 


0.86 ±0.22 ±0.12 



TABLE II: The values of the radius R in fm and the suppres- 
sion parameter A obtained by fitting the experimental corre- 
lation function with the model (2^ that takes into account 
the FSI effect in the resonance (/o +oo) approximation. The 
normalization = 1.03 in all cases. The values correspond to 
the third set of points from the right in Figure |7| so chosen as 
to strike a balance between statistics and purity. The results 
in the first and the second column respectively correspond 
to the single- and two-channel fits. The errors are, from left 
to right, statistical and systematic errors introduced by the 
uncertainty on the purity correction. The systematic errors 
from the model fits are very small in comparison and are not 
shown. 



eters used. Hence, the systematic errors are driven by 
the data and not theory. Figure |7| shows only a slight de- 
pendence of the radius parameter on the pair purity. On 
the other hand, A in panel (a) of Figure |H1 has a strong 
dependence on pair purity. Even though the purity cor- 
rection seems to improve the results, there is still a slight 
dependence remaining as shown in panel (b) of Figure 
IHl The value of A for the data with the highest purity, 
and therefore the cleanest signal, is consistent with unity. 
This is expected for a chaotic system with little contri- 
butions from decaying resonances. Plotting the radius as 
a function of the mean K^, as shown in Figure IHl shows 
a slight dependence of R with increasing Kt- However 
this could be a remaining artifact of the mean Kt de- 
pendence on pair purity, as mentioned earlier and shown 
in Figure 3. One has to look at several Kt bins for a 
specified pair purity to study a Kt dependence of the 
radius coming from real physics effects. This was not 
possible in this analysis due to the limited statistics. In 
order to strike a balance between statistics and purity, we 
averaged over the data from the coupled channel analy- 
sis corresponding the third set of points from the right 
in Figure [7{b), with a pair purity of « 82%, to obtain 
the values R = 4.09 ± 0.46(stat.) ± 0.31 (sys) fm and 
A = 0.92 ± 0.23(stai) ± 0.13(s?;s) at the mean transverse 
mass (tot) = 1-07 GeV. 

Figure ^Ul shows the tot dependence of R extracted 
from TTTT ||2^J, KgK'^, and proton- A correlations [s^. 
Considering the large mean transverse momentum of the 
pair, the value of R for K^ before taking into account the 
FSI in the K'^K^ system is larger than expected from the 



7p- 
6± 
5± 

3± 
2± 
1± 
7± 
6± 
5- 
4 

2^ 
1- 

8.5 



(a) 



Mass and coupling constants from Table I 
•a "b *c *d 



Ml 



(b) 



0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

Pair Purity 



FIG. 7: The extracted i? as a function of the pair purity (a) 
before correction for purity and (b) after correction for purity. 
The errors are only statistical. 



systematics followed by the rest of the data. However, af- 
ter taking into account the FSI effect the neutral kaons 
also seem to follow the tot scaling that hydrodynamics 
predicts [39l |. 
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Mass and coupling constants from Table I 
•a "b *c *d 



(b) 



5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

Pair Purity 



FIG. 8; The extracted A as a function of the pair purity (a) 
before correction for purity and (b) after correction for purity. 
The errors are only statistical. 




hb 0.8 0.85 0.9 0.95 1 1.05 1.1 

<K-,> (GeV/c) 

FIG. 9: The extracted i? as a function of the mean Kt of 
the pairs that go into the correlation function. The errors are 
only statistical. 



V. CONCLUSIONS 

We have presented the first measurement of neutral 
kaon correlations in heavy-ion collisions at RHIC. One 
has to consider the effects of FSI to obtain reasonable 
agreement between theory and data. The variations of 
the resonance parameters result in very small differences, 
which are well within our systematic errors. The effect 
of the pair purity on the correlation function has been 
studied extensively and is well understood. A Gaussian 



Kj, coupled-channel fit 



pA 
pA 
pA 
pA 



0.2 



0.4 



0.6 



0.8 



1.2 



1.4 



(GeV/c^) 



FIG. 10: i? as a function of rriT. Statistical and systematic er- 
rors are shown. The FSI uncertainty measured by the spread 
of the fit results in rows (a)-(d) of Table |^ is substantially 
smaller than the statistical error. 



fit to the correlation function does account very well for 
the C{Q) < 1 part of the data and gives a radius which 
is larger compared to the model fit results. 

The measured correlation radius is intermediate be- 
tween those obtained from two-pion and proton-lambda 
correlations in these collisions with the same conditions 
except for a different transverse mass, mx- The radii 
seem to follow a universal tot dependence in agreement 
with a universal collective flow predicted by hydrody- 
namics. The value of the parameter A, based on the high 
purity data, is consistent with unity and thus points to 
a chaotic kaon source. This is in correspondence with an 
indication of a dominantly chaotic pion source obtained 
from STAR measurement of three pion correlations ^3 ■ 

Our results represent an important first step towards 
a multi-dimensional analysis of neutral kaon correlations 
using the high statistics data from RHIC. In the future 
this analysis will allow to extract information about the 
freeze-out geometry, collective flow velocity, the evolution 
time and duration of particle emission. The latter is espe- 
cially interesting in the context of an increased emission 
duration expected if there is a first order phase transi- 
tion from a quark gluon plasma to a hadronic system. 
Recent pion interferometry measurements at RHIC how- 
ever point to a smaller evolution time and emission du- 
ration than expected from the usual hydrodynamic and 
transport models. This result may indicate an explosive 
character of the collision and is often considered as the 
interferometry puzzle. The fact that the Coulomb inter- 
action is absent in the dominant elastic transition and 
that the FSI effect can be handled with sufficient accu- 
racy makes neutral kaon interferometry a powerful tool 
which allows for an important cross-check of charged pion 
correlation measurements. Pion measurements are much 
more strongly affected by contributions from resonance 
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decays and final state interactions. 
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VII. APPENDIX 

The interaction of final state particles can proceed not 
only through the elastic transition ab — > ab but also 
through inelastic reactions of the type cd — > ab, where 
c and d are also final state particles of the production 
process. The FSI effect on particle correlations is known 
to be significant only for particles with a slow relative 
motion. Such particles continue to interact with each 
other after leaving the domain of particle production and 
their slow relative motion guarantees the possibility of 
the separation (factorization) of the amplitude of a slow 
FSI from the amplitude of a fast production process. For 
the relative motion of the particles involved in the FSI 
to be slow, the sums of the particle masses in the en- 
trance and exit channels should be close to each other 
[3^ . Thus, in our case, one should account for the effect 
of inelastic transition K^K~ — > K^K'^ in addition to 
the elastic transition K^K° -* K°K°. Instead of a sin- 
gle channel Scrodinger equation one should thus solve a 
two-channel one. In solving the standard scattering prob- 
lem, one should take into account that the FSI problem 
corresponds to the inverse direction of time. As a result, 
one has to make the substitution k* ~f ~k* and consider 
i^r'^if (= 1) as the entrance channel and K^K^{= 2) as 
the exit channel. Since the particles in both channels are 
members of the same isospin multiplets, one can assume 
that they are produced with about the same probability. 
Therefore the correlation function will be simply a sum of 
the average squares of the wave functions ^'"^^g, (t^) and 

\l/^^-^(r*) describing the elastic and inelastic transitions 
respectively. 

Assuming the s-wave dominance and r* outside the 
range of the strong interaction potential, one has '36f| : 

^\{n-ff{k*).[^^^p^, (12) 



where ^1 = m^o /2 and ^2 = /2 are the respective 
reduced masses in the two channels. p2 = ^2?'*, ''72 = 
{kla2)~^ and fc| = [2/^2 (fc*V (2^1 ) + Sm^o - 2mK+Y''^ 
is the momentum in the two-kaon rest frame. 02 = 
-(/i2e2)-i = -109.6 fm is the (negative) K+ K' Bohr 
radius, fc^{k*) is the s-wave transition amplitude re- 
normalized by Coulomb interaction in the K~^K~ chan- 
nel, G(p, 77) = ^ Ac (r/) . [Go (p, ?/) + iFo {p, 77)] is the combi- 
nation of the singular and regular s-wave Coulomb func- 
tions Go and Fq. Finally Ac{ri) — 2Trr] / [exp{27Tri) — 1] is 
the Coulomb penetration (Gamow) factor. 

The wave function of the elastic transition 1 — s- 1 is 
still given by Eq. in which k* = fcjf and the amplitude 
f ^ fl^ is now the element of a 2 x 2 matrix 

fc^(k-^ -ikr)"^ ■ (13) 

Here if is a symmetric matrix and kc is a diago- 
nal matrix in the channel representation: k],^ = fc*, 
k"^ = Ac{rj2)k2 — 2ih(T]2)l a2-, where the function h(r\) 
is expressed through the digamma function V'(^) = 
r'(z)/r(z) as hiji) = \^(i-r]) - ip{-iv) - \nr]'^]/2. As- 
suming that the isospin violation arises solely from the 
mass difference and Coulomb effects on the element fc^^, 
making it different from the momentum k* in the neu- 
tral kaon channel, one can express the K^^ matrix, in 
the channel representation through the inverse diagonal 
elements K^^ of the i^T-matrix in the representation of to- 
tal isospin I (the products of the corresponding Clebsch- 
Gordan coefficients being 1/2 or -1/2): 

{k-y^ = {k-^r = l[K-'-K-']. (14) 

The latter are assumed to be dominated by the reso- 
nances r = /o(980) and ao(980) for / = and 1, respec- 
tively, so: 

Ky^ ^{ml-s- ik'^-fDhr. (15) 

One should also take into account the correction 
AG^/f due to the deviation of the spherical waves from 
the true scattered waves in the inner region of the short- 
range potential, which is of comparable size to the effect 
of the second channel. This correction is also given in Ref. 
|36j| and is represented in a compact form in Eq. (125) of 
Ref. [33 ■ In our case, 

^^^^ = -i7^0/"W + l/c"W 

+ 21R{fl'fl'*)df] , (16) 

where dt^ = 2m{k~^Y^ /dk*^- at k* = 0, do coincides 
with the real part of the matrix of effective radii. 

One may see from Eqs. © and (|12|) that the usual 
resonance Breit-Wigner behavior settles only at small r* 
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when squares of the spherical waves \fl-' /r*]"^ dominate. 
At sufficiently large k* , one can neglect the Coulomb ef- 
fects and put fl^ ^ (/o + /i)/2, /^i = (/o - /i)/2, so 
that + P = l/oP + l/iP- The sum of the square 

terms then reduces to the incoherent Breit-Wigner con- 



tributions of /o and oq resonances. There can also be 
additional (not related to FSI) resonance contribution of 
the usual Briet-Wigner form due to direct /o(980) and 
ao(980) production. This contribution is assumed to be 
negligible as compared to the FSI effect. 
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